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1. Introduction. 

Exact solutions iiave a fascination that transcends any physical realisation. 
This is partly the reason for the use of toy, or cut down, models, the hope being that 
the generality and elegance of the solution would make up for any loss of realism. 
It is in this spirit that I have often resorted to evaluations in the Einstein Universe, 
which, although being static, is not a completely unrealistic cosmology, and it does 
form the instantaneous geometry of a Friedman-Robertson- Walker space-time. 

In some works, [1,2], I considered the Casimir effect on the Einstein Universe, or 
three-sphere, looking at the factored case, S"^/r, and, in particular, at fixed-point 
free, homogeneous actions. Because of the geometry of rotations, it is sufficient to 
evaluate for lens spaces, i.e. for F a cyclic group, Z^, and combine these to get the 
other cases. This has been done in [2] where attention was restricted to one-sided 
(homogeneous) lens spaces, the simplest kind. In the present work, after a question 
by Gregory Moore and David Kutasov, I wish to extend my treatment to two- 
sided lens spaces the reason being that these are also exactly computable, in one 
sense. The answers are no doubt academic and are presented in accordance with the 
remarks above. However, there is current interest in the possibility of a topologically 
non-trivial Universe and my results may be of relevance here. Relatedly, there has 
been recent discussion of the implications of the S^/F quotients for supersymmetry 
in the context of the AdS/CFT correspondance and the brane world scenario [3]. 

Because the physical motivation is meant only to lead quickly to the mathe- 
matical calculation, I concentrate on the total vacuum energy, as opposed to the 
energy density. In the next section some basic facts are outlined that lead as quickly 
as possible to the technicalities, which, I admit, are my main interest in this work, 
which should be considered as a continuation of [2]. 

One reason for the attraction of the three-sphere is that it is a group manifold. 
Because of the isomorphism S"^ ~ SU(2), it is handy to implement the action of 
F on the sphere as the combination of left and right group actions on SU(2) i.e. 
roughly, F ~ Fj^ x F^. A general point on the sphere S^ is denoted by q, or, as an 
element of SU(2), again by q (or sometimes by g) and the general action g — > 75 is 
realised as g — > IlIIr- 



2. Modes. Casimir energy. 

I intend to consider (massless) fields of spin 0, 1/2 and 1. The spin zero 
expression for the Casimir energy has been given in [1], and repeated in [2], for 
double-sided actions. For one-sided actions, the corresponding formulae for spins 
1/2 and 1 have been detailed in [4,5] and [6]. For completeness I wish to derive 
the double-sided forms. The answers could be be deduced from the off-diagonal 
^-function given in [7] but it is preferable to start again. Unfortunately, a certain 
amount of dynamic and, worse, kinematic recapitulation is required to reach the 
needed forms. 

The equations of motion that concern me here are the Klein-Gordon, Dirac and 
Maxwell equations. For mathematical convenience, I consider conformally coupled 
scalars so that the eigenvalues of the relevant Helmholtz equation are squares of 
integers. 

In these cases there are no field theory divergences to worry about, in the ab- 
sence of fixed points, and a reasonably unambiguous formula for the Casimir energy 
can be given as a regularised sum of the eigenvalues of the relevant Hamiltonian 
operator which, for scalars, is the square-root of the Helmholtz operator. The 
spectrum has positive and negative parts. In order to be more specific, some mode 
properties are required and I now briefiy repeat some standard kinematics, [8,7], 
firstly on the full sphere. 

For spinor fields (i.e. any field with spin greater than zero), the Cartan moving 
frame method projects everything onto one of the sets of Killing fields on S^. It is 
convenient to use the left-invariant set and we then say we are using 'right spinors'. 
These are scalars under left actions and, further, the standard results of angular 
momentum theory can be used with no changes in standard conventions. 

For spin-zero, the angular momentum representation matrices, l)^^{g), L = 
0, 1/2, 1, . . ., provide a complete set of modes on the full three-sphere. 

For higher spins, j (= 1/2, 1, . . .), spin-orbit coupling gives the modes as the 
spinor hyper spherical harmonics, 

Yifj{q) = {m,q\JMLjN) 

(2L + l)(2J + l)V/'^^^,^^^^j L M\ (1) 



27r2a3 J '' '^\mN' J 

The equations of motion imply the massless polarisation conditions 

J = L ± j for L > j 
J = L + j for L < J , 



and then it can be shown by simple spin-orbit diagonahsation that the eigenvalues 
of the Hamiltonian operators (the square-root of the Helmholtz operator, the Dirac 
operator and the Maxwell curl) are, 

H I JMLjN) = El J I JMLjN) , (2) 

with 

The energy spectrum on the full sphere is thus 

E+ = -{j + L), 1=1,2,... 

a (3) 

EZ = -(j-L), L = 2j + l,2j + 2,.... 
^ a 

1 note that E~^ is positive and E~ is negative, and that the spectrum on the full 
sphere is symmetrical about zero. 

As is well known, factoring the sphere to S^/F modifies the spectrum by elim- 
inating certain modes and the spectrum is a subset of the full sphere one. I denote 
the resulting degeneracies of the eigenvalues (3) by d~^{L) and d~{L) and define the 
handy positive and negative Hamiltonian ^-functions. 
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(4) 



in terms of which the regularised value of the Casimir energy is 

E, = {-if'\h{j) (c;(-i) +c-(-i)) , (5) 

where h{j) is a spin degeneracy factor, h{0) = 1 and h{j) = 2, j > 0. Corners 
have been cut here in anticipation of the fact that there is no divergence in the 
continuation of the (^-functions at —1. 

The problem now falls onto the evaluation of the degeneracies which can 
be deduced in the following rather oblique way. Rather than analyse the projected 
modes, it is clearer to consider the covariant heat-kernel for, say the relevant second 
order spin-j propagation equation of which (1) are eigenmodes and the eigenvalues 



are the squares of Ej^ j, (2). The associated degeneracies wiU be denoted by d{Lj J) 
and are related to those in (4) by, 

d+{L) = d{L,L + j) , d-{L) = d{L,L~j) . (6) 

The spinor heat-kernel on S^/F, Kp(g, q'), is given, as usual, in terms of that, 
K.{q,q'), on S"^ by a pre-image sum. However, because right spinors have a non- 
trivial behaviour under right actions, it is necessary to provide compensating right 
'gauge' rotations at the pre-image points to bring the beine into uniform alignment 
so that the pre-image sum makes sense. 

One is therefore interested in the traced heat-kernel on the factored space, 

^r= E / c^qIV[K(q,7^q7^)P^(7^1)], (7) 

where K is constructed from the eigenmodes (1). In order to use the properties of 
these, the integral in (7) is changed to one over the full sphere by using the invariance 
of the volume, dq, under symmetry actions and also the equation, following from 
the definition of right spinors, [8], 

K{^qV.Cq'v)=1^'iv-')^iq.q')1^'{v)- (8) 

Then 

^r = ^ E / dqTr[K{q,^^q^^)V\^-')]. (9) 

The switching equation, (8), shows that it is sufficient to know K((7) = K.{g, id) 
because K.{g',g) = 'K{g~^g',id) (note the order). 

The construction of the heat-kernel, K, in terms of the modes, Y, gives us the 
expression for the degeneracies {\M.\ = 27r^a^), 

d{L,J) = 

Group products and the orthogonality of the P's easily yield 

/ dqV^^,^ij^q-'^,q) = ^X^^^(7j*^(7«) 



and then standard angular momentum manipulations in the end give simply, 

d(L,J) = ^^ J2 X^''\il)x^'\ir) (10) 

in terms of the SU(2) characters, 

, J X sin L9^ 

smfc/^ 

where 6 is the 'radial' angular coodinate labelling the group element, 7. 

For spins greater than zero the spectrum is no longer symmetrical about zero, 
as is well known. 

3. The Casimir energy formulae. 

The equations (5), (4), (6) and (10) can now be put together to give explicit 
summations for the Casimir energy. Thus 

(-1)^^\/.N i:..„ „s V^ Ty^ 1 sinn^^sin(n + 2j)^^ 



J 4r ^-^'s^-i ^-^ A^tn+jY sm6'rsm6'o 

1 sin n6i^ sin(n — '2.j)6^ 



y 

^-^ (n — j)^ sin^r sin^R 

n=2j + l V J/ L H 

where I have set n = L and 9j^ = 9^ , 6j^ = 9^ . 

It must be admitted that this equation is meant to apply only for j = 0, 1/2 
and 1. The further evaluation is possible for all j but the field theoretic significance 
is unclear. 

The only n-sum in the group sum to diverge as s -^ — 1 is that for 7 = id, 
i^R = ^L = 0)- So this term is separated and the limit taken in the rest. 

The identity contribution is 



Ef = ^ , \ h(j) hm a' 

J 4r ^ s^-i 



y> n{n + 2j) y^ n{n - 2j) 
^ (?i+ ?')^ ^ (n- jY 



which can be evaluated as Hurwitz (^-functions giving, expanding the Bernoulli 
polynomials, (see [7]), 



Specific values are 



rpid rpi' 



17 



/2 960a|r| 



Ei^ 



11 

120a|r| ' 



(12) 



which, apart from the volume factor, l/|r|, are the known full-sphere values. 

The harder part of the calculation is the evaluation of the non-identity contri- 
bution, 



E E(^+-^' 



..sinnOj^ sin(n + 2j)9j^ 



7i,7«'-n=l 



sin 6*^ sin 6*^ 



n=2j + l 



., sin 716*^ sin(n — 2j)9j^ 



sin 9]^ sin 6*^ 



(13) 



It turns out that the integer and half-odd-integer j cases differ and so I treat 
them separately. For integer j, a shift of summation variable yields the more sym- 
metrical form, 



E' 



1 



E 



sin 6*1^ sin 6*0 ^^-^ 



n 



sin(n — j)9j^ sin(n + J)(^r + sin(n + j)9j^ sin(n — j)9j^ 



and simple trigonometry produces. 



-— > n f cos nd cos ?'a — cos na cos ?'/9) . 

sin 0r sin 6*0 ^^ \ "^ -^ -^^J 

11 '7fl ^ ^ n=j 



where 



a = 9j, + 9L, P = 9j,-9l. 
The sum over n can be done, 



-^E' 



a,P 



cosec ^/3/2 cos ja — cosec '^a/2 cos j 13 
cos a — cos (3 



While this expression can be taken further for any j, 1 now set just j — and 
j — 1 for simplicity, when again simple algebra gives 

cosec ^/3/2 — cosec ^q;/2 



--E' 

a, 13 '- 

for J = and 



cos a — cos 13 



-^ cosec ^/3/2 cosec ^a/2 , (14) 



a, 13 



2^ 

ot,l3 



cosec ^/3/2 cos a — cosec '^a/2 cos /3 
cos a — cos 13 



-^ (cosec ^a/2 + cosec ^/?/2 cosec ^a/2 cosec ^/3/2) 



(15) 



a, (3 



for j = 1. 

For spin-half I return to (13) and rewrite the sum as one over odd integers, 
n = 2n + 1. Trigonometry and use of the a and /3 now give 



E' 



1 



cos a — cos (3 
and the n summation is done using 



V^ n (cosn/3/2 cosq;/2 — cosnQ;/2cos/3/2) 



oo ^ oo 

2_] smn6 = -cosec6', \J 



n=l,3 



n COS n9 = — cot 9 cosec 9 . 
2 



which yields, 



n=l,3 



- y^ cot a/2 cosec a/2 cot 13/2 cosec /9/2 . 



(16) 



4. The Casimir energy calculated on lens spaces. 



Reinstating the factors, the expressions for the three Casimir energies are 



-^1/2 

El 



1 



a\T\ 
1 

a\r\ 

1 

airi 



1 1 



240 16 



y cosec ^/3/2 cosec ^a/2 



a,/3 



1 — > cot q;/2 cosec q;/2 cot /3/2 cosec /3/2 

960 s^^ ' / y/ y/ 

11 1 



«,/3 



120 4 



+ T^^ (cosec a/2 + cosec /3/2 — -cosec ^q;/2 cosec ^/3/2) 



«,/3 



(17) 

The expression for Eq was given in [1]. For one-sided actions (i.e. a = ±/?), 
the expressions reduce to the ones in [2,6]. 

For the lens space, L{q; l^, I2), the angles a and (3 take the values 



27ipiy^ 2Tcpiy2 

a= , [3 = 



Q 



Q 



(18) 



where p, = 1,...,^ — 1, labels 7. Ui and z/2 are integers coprime to g, with li and I2 
their mod q inverses. The simple, 'one-sided' lens space, L{q; 1, 1), corresponds to 
setting U2 = i^i — ly' — I, say, so that 9j^ — and 6*^ — 2Tvp/q. These were discussed 



in [2] . I now examine the general case and begin with the integer spin forms which 
involve only powers of cosecants. The single sums of cosec^ are classic, 

S{q; ^) = -y^ cosec = , (19) 

and independent of z/, being sums over the roots of unity ordered in different ways 
as in the stellated polygon. 

I require the harder, split summations, 

6(0; z/i , z/n) = - > cosec cosec (20) 

^Pl Q q 

where u^ and 1/2 are fixed integers from 1 to g — 1. For these particular forms 
I can fortunately make use of the computations of Zagier [9], who has defined a 
generalised Dedekind sum (I change his notation slightly) 

r(,; .„.„...,. J = (-l)V2iX: cot ^^... cot ^^. (21) 

p=i 

Here g is a positive integer, z/j^, . . . , z/„ are integers prime to q and n is even. I note 
the important fact that one of the sets of roots of unity can be reordered to the 
canonical one while maintaining the relation between all the sets so keeping the sum 
the same. Formally this is done by multiplying all the u^ by the mod q inverse, Z^, 
say. 

I need the four-dimensional case, r(g; z/^^, z/^^, z/2, z/2), and, as said, one can put 
z/]^ = 1 without loss of generality. So I look at r{g; 1, 1, z/, u) and note the simple 
relation that takes me back to the cosecant sums, 

Siq; z/, 1) = r{q; 1, 1, z/, z/) + 2S{q; 1) - ^^ 
= T{q; 1, 1, z/, z/ + ^ ^p ^ , 

S{q] 1) being given by (19). 

In his Table 3, Zagier has listed some four-dimensional examples, of which 
(i(q'; 1, 1, 3, 3) and d{q; 1,1, 4, 4) are most relevant for us. The first quantity has 
also been evaluated by Harvey et al, [10], in its cosec version. Employing Zagier's 
expression, one finds, 

Q( Q 1^ g" + 210q^ ± 80q - 291 

S{q;?,,l) = — , Q = ±l(mod3) 



in agreement with [10]. Note that use of the cosec form has resulted in a more 
symmetrical result. 

My intention here is to take the computations further for r(q'; 1, 1, z/, z/) with 
larger v. 

5. Evaluation of the Dedekind sums. 

The essential calculational tool is the (n + l)-term reciprocity law obeyed by 
the Dedekind sums (21) (I state it for the 4-dimensional case), 

4 

^t(z/j;z/o,...,z^j,...,z^4) = </'4('^o.---.'^4), (23) 

where z/q, . . . , 2/4 are pairwise coprime positive integers, the hat signifies omission of 
the indicated term and (f)^ is a function whose form can be calculated. 

Zagier derives the relation (23) using contour integration which indicates why 
the construction of r in terms of products of cotangents is so convenient as the pole 
defining function is another cotangent, giving a symmetrical structure. This leads 
to the expression for (j)^, 



(t^n{^0^---^^r> 



^ ^n{^0^---^^n) 



^0---^n 



where the polynomial, £^(z/q, . . . , z/^) is the coefficient of t" in the power series 
expansion of 

-'■-'- tanhz/-t 

The individual expansions of the coths give an explicit multiple sum expression for 
this quantity (Berndt, [11]) which is better expressed as a generalised Bernoulli 
polynomial but I do not need this cosmetically compact form here. 

The relation (23) cannot be applied immediately to the case here because the 
members of the set, (^,1,1,//,//), are not mutually coprime. However z/ is prime 
to q so, following Zagier, I make use of periodicity to replace T{q;l,l,u,u) by 
— r(q'; 1, 1, z/, g — u) so that the set (5,1,1,//,^ — z/)) is coprime and (23) can be 
applied. In writing (23) out I immediately use the fact that r(l; z/]^, z/27 ^^3? ^^4) is 
always zero which produces a recursion formula, as is now shown. The reciprocal 
relation reduces directly to 

T(g; 1, 1, z/, g - z/) + t{u] 1, l,q,q-iy)+ r{q - z/; 1, 1, z/, q) = (p^{q, 1, 1, z/, g - z/) (24) 



For short, set 

so the first term on the left in (24) is — r and the third term is t ^^ as can be shown 

using mod q — u periodicity. The second term is to be taken to the right-hand side 

and together with (j)^ will form the driving term of the recursion which, as is seen, 

successively reduces q by z/. The driving term has to be explicit, (f)^ is given by 

Zagier and equals, in this case, 

, , , , , , 3z/4 - 6z/3g + 9z/2g2 + 20z/2 - z/q3 _ 20z/g + 3^^ + 20^^ + 3 
<^4 W, 1, 1, z/, Q - z/) = 1 — — ^ 

(25) 

A little more argument is needed for r(z/; 1, 1, g, g — z/). It is necessary to break 

up the solution according to the residue classes mod z/ of q. Thus, if g = nz/ + x, 

r(z/; 1, 1, g, g — z/) = r(z/; 1, 1, x, x) and this 'constant' has to be calculated separately 

for each x, (1 < a; < z/), for the given z/. 

As an introduction to the general method, I check the case mentioned earlier, 
z/ = 3. There are two ranges for q, i.e. q = ±1 mod 3. One therefore requires 
r(3; 1, 1, 1, 1) and r(3; 1, 1, —1, —1), = r(3; 1, 1, 1, 1) also. This quantity follows from 
the classic summation, 

r 1 1 1 1^ IV- ,4^P (g-l)(9-2)(g^ + 3Q-13) 
r(Q;l,l,l,l) = -^cot — = — , (26) 

as 2/27. Hence, evaluating the function ^4, (25), the recursion reads, 

142 _ 142 q^ q 17 

^-^ ^ 135^ ~ ^'^-^ ^ 135(q-3) ^ 45 ~ 15 ~ 45 ' 
in both cases. This can be iterated down to r^^ = for g = 1 mod 3 and to T2 = 
r(2; 1, 1, 3, 3) = r(2; 1, 1, 1, 1) = 0, also, for q = -1 mod 3. Thus for q = 3[q/3] + 1, 

- _ill V (^(^^±11! _ 3fc + l _ 17\ 142 
^'^ ~ 135q ^ ^ V 45 15 45/ ^ 135 

{q - l){q^ + q^ - 59q + 426) 



405g 

and for q = 3 [5/3] + 2, 



142 'y^Y(3A; + 2)2 3^ ^ 2 17\ 142 
^'^ ~ "135^ ^ ^ V 45 15~ " 45 J ^ 270 ^27) 

{q - 2){q^ + 2q^ - 56q + 213) 



405q 
which are the forms listed in Zagier, [9]. 
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6. Automatic evaluation. 

The above 'hand' calculation is typical. I wish to see how far it can be auto- 
mated for any, but fixed, p. 

The form of the function (f)^ allows one to extract the combination, 

3z/4 + 20Z/2 + 3 
"^ = "^ + ^^q ' ^''^ 

and to write the general recursion as, 

^ , 3q^-3z/q + 6z/^ + 20 
(Jq - (Tq-^ = r{iy; 1, l,x,x)-l + — . 

For q = \qlv\v + x, for given x, this is taken down to o^ to give ^ 

[q/^] r 



'^q =^2 






3(z//c + x)^ - 3z/(z/A; + x) + 6z/^ + 20 



r(z/;l,l,a;,a:)- 1 + 

45z/ 



(T^ + (r(z/; 1, 1, x, x) - \) — H j^^ 



(29) 



V 45z/^ 

Everything is therefore explicit except the starting value, a^, and the 'constant', 
r(z/; 1, 1, x, x) which have to be found separately. The tables in Zagier could be 
consulted after suitable transformations using the basic properties of the r's to 
reduce the arguments to those appearing in the table but I here continue with the 
recursion approach for a self-contained treatment. 
Some streamlining of notation is indicated. I set, 

Tq = T(g;l,l,'^, z^) = T{q,iy), 
and rewrite the iteration result (29) in terms of the r's, 

T{q, u) = r(x, u) + ^i^ r(z/, x) + F{q, u, x) , (30) 

V 

where, 

1 



F(a, 6, c) = ^ ,; ( a(36^ + 206^ + 3) - ac^ - hac^{}? -96 + 4)- 



c(36^ + 56^(4 - a^) + 45ba^ - a^ - 20a^ + 3) j 



(31) 



X does not alter during the iteration. 

11 



I recall that q is arbitrary, z/ is a fixed number, prime to q and 

X = qvaodp (32) 

is to be considered a chosen number, also prime to q. 

In the first term on the right-hand side of (30), u can be reduced mod a;, 
defining, 

X = V mod X , (33) 

[x' < X <v < q). Hence follows the fundamental relation, 

T{q, v) = r(x, x') H riu, x) + F{q, z/, x) . (34) 

I introduce the series of Euclidean relations, 

a^i =Xi_2modXi_i, x^_^<x^, (35) 

of which (32) and (33) are the first two examples, with Xq = q, x^ = v, Xi = x and 
0:3 = x' . The recursion (34) then becomes, with a more systematic labelling, 

t(xj, x^^^ = T{Xi_^_ij 3;j_|_2) + r(a:j_|_2, a^i+3) + -r (a^^i, a^j+i, Xi_^2) ^ (36) 

which begins with z = giving (34). The task, to repeat, is to find the 'constants' 
t{xi,X2) and t{x2jX^) by iterating (36). 

The algorithm (35) ceases when one of the Xj reaches 1. Assume that this 
happens when i = t {t is usually small). The iteration, (36), will also come to an 
end when i + 2 = t because r is known whenever one of the arguments is 1 , 

{b-l){b-2){b' + 3b-13) 

^^^' '^ = 456 (37) 

T{l,d) =0. 

This means that t{x^_2,x^_i) is known from the recursion which can be rolled 
up from the bottom by successive back substitution. I find, for example, the first 
non-trivial value, 

, ,, a-1 (6-l)(6-2)(62 + 36-13) ^, , , 

_ (a - 1) (a3 + a2 + a{b^ - 15b^ - 5) + 36^ + 206^ +3) ^ ^ 

~ 45a62 ' 

12 



where (a, 6, 1) are the last three numbers in the Euchdean algorithm, (a:j_2, a^t_i, a^t)- 
In particular, if a = 6 + 1, 

^('+'^'^- 45(6+1) - 

which also follows simply from the reduction, 

r{b + 1; 1, 1, b, b) = rib + 1; -6, -6, 6, 6) 

= r(6 + 1; -1, -1, 1, 1) = r(6 + 1; 1, 1, 1, 1) , 
and (26). 

Likewise, but more generally, one finds (see also the Appendix), 

r(„6+l;l,l.n,n) = J 45(iTTr ' ''" 

Similar reductions have to be performed in order to reach one of the given 

forms in Zagier's table. My formula has the manipulations built in. 

The systematic elimination is therefore performed by inverse recursion and a 

more compressed notation is useful. I define 





«fc 


_ ^t+i-k ~ 

^t+2 


-k 


k 








^k 


= '''i^t+l-kJ 


^t+2- 


k) 








fk 


= F{xt^^_k 


7 ^t + 2- 


-k^^t+3- 


-.), 




so that the two-term 


recursion (36) becomes 










T-fc = afcTfc.i + Tfc. 


2 + fk 


, k = 


3.. 


.t. 


with the starting values. 












^1 -- 


= 












^2 -- 


{h- 


l){b-2){9 


+ 36- 


-13) 

7 


6 = 


--Xt-l 



(40) 



To make everything plain, my two step procedure is to compute the quantity 
I want T{q, u) = t{xq, Xi) from (36) at z = by finding the constants t{xi, X2) and 
r(x2, X3) from the recursion (40), which has to be taken to k = t — 1 and k = t. The 
reason is that this second stage is purely numerical. The variable q is arbitrary, 1/ 
is a constant and a; is a number chosen in the range 1 < x < u. (It may be possible 
to leave x variable, but I prefer to set it by hand.) 

The first act is the standard computation of the Euclidean algorithm (35) for 
given Xi = 1/ and X2 = x. It is then straightforward to implement the recursion to 
produce the polynomials in q in one go. I coded the entire process quite compactly 
in DERIVE. 

13 



7. Some results. 



The general expression for the cosecant sum, q S{q; z/, 1), is a polynomial which 
I write in the following way, 



q S{q- V, 1) = A(/ + Bq^ +Cq + D) 



(41) 



where, introducing x = q mod z/, the coefficients A, B, C, D depend on u and x. 
A, B, D are unchanged under x ^ u — x, while C changes to sign. So x can be 
restricted to 1 < a; < z//2, and, of course, x and u are coprime. 

When X = 1 or 2, there are relations among the coefficients because the cotan- 
gent polynomials vanish when g = 1 or 2, respectively. 

The factor of A has been extracted because it is known from the obvious g — ^ oo 
limit of (20), 3 



A 



, Cfl(4) 
'7r4z/2 



1 



45z/2 



(42) 



(The factor of 2 is a consequence of the periodicity of the summand.) Hence only 
B, C and D need be displayed which I do in the following table for a few low values 
of u and x, 



V 


X 


B 


C 


D 


V 


X 


B 


C 


D 


3 


1 


210 


80 


-291 


9 


1 


7770 


12320 


-20091 


4 


1 


490 


360 


-851 


9 


2 


3450 


4960 


-20091 


5 


1 


994 


1008 


-2003 


9 


4 


3450 


-640 


-20091 


5 


2 


706 


144 


-2003 


10 


1 


11494 


19008 


-30503 


6 


1 


1830 


2240 


-4071 


10 


3 


4006 


3456 


-30503 


7 


1 


3130 


4320 


-7451 


11 


1 


16450 


28080 


-44531 


7 


2 


1690 


1440 


-7451 


11 


2 


6370 


12240 


-44531 


7 


3 


1690 





-7451 


11 


3 


4930 


6480 


-44531 


8 


1 


5050 


7560 


-12611 


11 


4 


4930 


3600 


-44531 


8 


3 


2170 


1080 


-12611 


11 


5 


6370 


-2160 


-44531 



The lowest case, z/ = 2, for which B = 70, C = and D = 71 can be obtained 
directly using the relation cosec^26' = (cosec^6' + sec^6')/4 and the classic, single 
summations of cosec^, cosec^ and tan^. 



■^ This can be applied also to the general sum, (21). 
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8. Some Casimir values. 

In this paper I consider only the integer spins, and 1. The basic expressions 
for the Casimir energies are given in (17). 

In terms of the sums discussed above, the spin formula on the lens space 
L{q; 1, 1), where / is the mod q inverse of u, is. 



Eoiq, ^) = - 



- T^S{q]P,l] 



240q 16 



and for spin-1. 



EA<i.^) = Ya 



Wq^S{q-l)-\s[q-.A) 



(43) 



(44) 



As a sample, I plot in fig.l the spin-0 energy as a function of the order, q, for 
a fixed twisting, z/ = 5, 



X 



+ 1 « 2 o 3 X 4 




figl. Casimir energy E for conformal scalars on lens spaces of order q and 
twisting p — b. x — q mod 5 
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Fig. 2 plots, typically, the Casimir energy against the twisting v for a fixed order 
q = 29. It exhibits the symmetry about the value u = q/2 and a series of minima 
the significance of which I do not know. The spin-one graph is similar except it 
goes positive for certain values of u. 




twist 

fig2. Casimir energy E for conformal scalars on a lens space of order 29 for 
twistings z/ = 1 to 28. 
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9. Conclusion and discussion. 

The values for the Casimir energies presented here, and in our earlier works, are 
doubtless unpractical but their exact determination is not without methodological 
value. 

Harvey et al, [10], present an algorithm for the computation of the cosecant 
sums, which can also be automated. In content, it is, of course, the same as that 
given here. However I feel the present approach via the cotangent sums is preferable 
as being more explicit and systematic. 

There exist many generalisations of the original Dedekind sum. Some of the 
older ones are covered in the book by Rademacher and Grosswald, [12]. A more re- 
cent sum, extending that of Zagier, has been proposed by Beck, [13], which involves 
replacing each cotangent in (21) by a higher derivative. Beck employs contours 
to determine the reciprocity relation which could be used to compute the cosecant 
sums directly by simply restricting to first derivatives. I have not checked this. The 
article by Berndt and Yeap, [14], contains further references. 

Unfortunately I could not find the spin-half summations, (16), so their recip- 
rocal relation will have to be sought using contours, ab initio. 

Also, as a challenge, if one could evaluate the even higher expressions, 

rfg; 1, 1, z/, z/, . . ., u, u) = - > cot — cot ...cot , (45) 

n •^— — ^ n n n 



the cosecant sums, 

- > cosec — cosec . . . cosec , (z/ , . . . , (U ) C (z/, . . . , ^u) , (45) 

Q^ q q q 

would follow on binomial expansion. Such sums would be needed when evaluating 
Casimir energies on the higher spheres. 

For the case here, a more careful independent analysis of the asymptotic be- 
haviour of the cosecant, or cotangent sums as g -^ oo would yield general informa- 
tion about the coefficients in (41). 
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Appendix. 

I make some remarks on the evaluation of tiie Dedekind sums not necessary 
for tiie immediate prosecution of the calculation of the Casimir energies but which 
extend the results a little. 

In the preceding calculation the Euclidean algorithm was taken down to its 
x^ = 1 limit, which always works. However it is possible to stop if the penultimate 
value, x^_i, equals 2 because r is also known when one argument is 2, 

T{xt_^,Xt) =r(2, 1) = 

(a-l)(a3 + a2-49a + 131) 



r(a,2) 



180a 



The second expression checks against the result of the earlier recursion, and is in 
Zagier's list, but the point is that it follows directly from the standard sums of 
(even) powers of cotangents (or of cosecants) after using the simple trigonometric 
formula, 

2 cot e cot 29 = cot^ e-1. (47) 

For convenience, I denote the last four entries in the Euclidean algorithm, 
(a;i_3, Xf_2, a^t_i, x^), by {d, a, b, c) with 6 = 2, c = 1. The algorithm implies that a 
is odd, a = 2n + 1, and also that d = ma + b, (n, n G Z). 

The first non-trivial result of the backwards recursion is 

T{d, a) = T{a, b) + F(d, a, b) , 

a 

or, putting in the values, 

r(m(2n + l) +2,2n + l) 

TTl 

' m^(4n^ + 4n + 1) + 8m^(2n + 1) 



45(m(2n + l)+2) 

+ m(4n^ + 8n^ - 24n^ - 28n + 4) + 16n^ + 24n^ + 2n - 3 

Setting n = 1,2,... reproduces a subset of our previous expressions. For example 
n = 1 gives (27). The advantage of this form, and its companion, (39), is the explicit 
dependence on n. 

Incidentally, the identity (47) can be used to obtain the specific value of the 
original Dedekind sum, s{2,q) = {q ~ l){q — 5)/24q. More generally it allows 
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one to obtain explicit forms for the sums, T{q; 1, . . . , 1, 2, . . . , 2), as combinatorial 
combinations of generalised Bernoulli polynomials in q yielding, for example, 

T{q; 1, 1, 1, 2, 2, 2) = -^^{q -l){q- 5)(2q4 + Uq^ - 43q'' - 318q + 995) 

T{q; 1, 1, 1, 1, 2, 2) = -j^{q -l){q- 2){q - 4){2q^ + 7q^ -7q- 139) 
t(q;1,1,1,1,1,2,2,2,2,2) = 



;^— — (q -l){q- 5){2q^ + 12q'^ - 103/ - 678g^ 



2993760Q 

+ 20139^* + 15468^^ - 18017^2 - 185442^ + 344425) . 

Curiously, the symmetrical quantities seem to vanish when g = 5 on the spheres S"^, 
S^\ Si9 ...but not on S^ S^^ .... 

The identity also permits a sum such as r(g; 1, 2, z/, . . . , ^u) to be replaced by the 
'simpler' ones, r(g; 1, 1, z/, . . . , ^u) and r(g; z/, . . . , (u). More pertinant for this paper, 
it means that r(g; 1, 2, z/, z/) can be reduced to r((3'; 1, 1, z/, z/) {q odd), computed 
previously, and to r(g; z/, u) = T{q] 1, 1), a standard sum, 

r(g; 1, 2, z/, z/) = lr(g, z/) - l(g - l)(g - 2) . 
2 D 

The table in Zagier enables one to check this relation, which I would like to regard 
as extending the scope of my earlier results. 

The consequences of more complicated trigonometric identities are discussed by 
Zagier, [9]. He shows, for example, that r(g; 1, 1, 2, 4) = r(g; 1, 2, 2, 2) by an identity 
following from (47). Actually, repeated direct application of just (47) reduces each 
of these quantitites ultimately to standard summations, so rendering Zagier's earlier 
computation via reciprocity unnecessary. 
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